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Abstract 

An elegant formula for coordinates of carbon atoms in a unit cell of a single-walled nanotube 
(SWNT) is presented and a new molecular motor of double-walled carbon nanotube whose inner 
tube is a long (8,4) SWNT and outer tube a short (14,8) SWNT is constructed. The interaction 
between inner an outer tubes is analytically derived by summing the Lennard-Jones potentials 
between atoms in inner and outer tubes. It is proved that the molecular motor in a thermal bath 
exhibits a directional motion with the temperature variation of the bath. 
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It is well known that a pollen in water exhibits Brownian motion. The forces on the pollen 
stem from two components fl]: A fluctuating force that averages to zero over time, and a 
viscous force that slows the motion. These two kinds of forces are related by temperature 
so the fluctuation is often called thermal noise. The second law of thermodynamics suggests 
that biased Brownian motion requires two conditions Q]: (1) breaking of thermal equilibrium 
and (2) breaking of spatial inversion symmetry. In order to illustrate the second law of 
thermodynamics can not be violated, Feynman devised an imaginary ratchet system with a 
pawl in his famous lecture Q]. 

To explain the matter transport in biological systems, the concept of molecular motor is 
introduced 14, nf. From the statistical viewpoint, many models |2j of molecular motor were 

no n n 

put forward, such as on-off ratchets [y, 12| , fluctuating potential ratchets |8| , fluctuating 



force ratchets 0, [l(3] , temperature ratchets [ill 



12l [l3j and so on. All those models satisfy 



the two conditions required by the second law of thermodynamics. 

With the development of nanotechnology, especially the discovery of carbon nanotubes 
Q|, people are putting their dream of manufacturing nanodevice [3] into practice. Here a 
natural question arises: Can we construct molecular motors from carbon nanotubes? In this 
Letter, we will prove that it is possible to construct a molecular motor from a double-walled 
carbon nanotube (DWNT). The molecular motor in a thermal bath exhibits a directional 
motion with the temperature variation of the bath. 

A DWNT consists of two single-walled carbon nanotubes (SWNT's) with a common axis. 
The layer distance between the two tubes is about 3.4 A. A SWNT without two end caps 
can be constructed by wrapping up a single sheet of graphite such that two equivalent sites 
of hexagonal lattice coincide |l6[. To describe the SWNT, some characteristic vectors need 
to be introduced. As shown in Fig^ the chiral vector C^, which defines the relative location 
of two carbon atoms, is specified by a pair of integers (ni,n 2 ) which is called the index of the 
SWNT. We have Ch = Wiai + 722^2 with ai and a2 being two unit vectors of graphite. The 
translational vector T corresponds to the first lattice point of 2D graphitic sheet through 
which the line normal to the chiral vector Ch passes. The unit cell of the SWNT is the 
rectangle defined by vectors and T, while vectors &i and a 2 define the area of the unit 
cell of 2D graphite. The number N of hexagons per unit cell of the SWNT is obtained as a 
function of n\ and n 2 as N = 2{n\ J \- n\ + nin 2 ) / where d^ is the greatest common divisor 
of (2n2 + rix) and (2ni + n%). There are 2N carbon atoms in each unit cell of the SWNT 
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because each hexagon contains two atoms. To denote the 2N atoms, we use a symmetry 
vector R to generate coordinates of carbon atoms in the nanotube; it is defined as the site 
vector having the smallest component in the direction of C^. From a geometric standpoint, 
vector R consists of a rotation around the nanotube axis by an angle — 2tt/N combined 
with a translation r in the direction of T; therefore, R can be denoted by R = (^(t). Using 
the symmetry vector R, we can divide the 2N carbon atoms in the unit cell of SWNT into 



two classes: one includes N atoms whose site vectors satisfy 
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A, = ZR- [ZR- T/T 2 ]T (Z = 0, 1, 2, • • • , N - 1), (1) 
another includes the remainder N atoms whose site vectors satisfy 

* = m+Bo -[fiK^]T 

_ [ (iR + B ).C > ]Ct (( = 0ili2 _... iAr _ 1)i (2) 

where B represents one of the nearest-neighboring atoms of A . In and only in the above 
two equations, [• ■ ■] denotes the Gaussian function, e.g., [5.3] = 5. 

We construct the molecular motor of double-walled carbon nanotube as shown in Fig|2j 
The inner tube's index is (8, 4) with a length long enough to be regarded as infinite. The 
outer tube is set to be a (14, 8) tube with just a single unit cell |l8| . Obviously, they are 
both chiral nanotubes and their layer distance is about 3.4 A. If we prohibit the motion 
of the outer tube in the direction of nanotube axis, it will be proved that this system in a 
thermal bath exhibits a directional rotation when the temperature of the system varies with 
time. Thus it could serve as a thermal ratchet. 

To see this, we first select an orthogonal coordinate system shown in FigEl whose z-axis 
is the tube axis and x-axis passes through carbon atom Aq of the inner tube. We fix the 
inner tube and forbid the ^-directional motion of the outer tube, that is, carbon atom A of 
outer tube is always in the xy-p\ane. We denote the angle rotated by the outer tube around 
the inner tube to be 9. 

We take the interaction between atoms in the outer tube and the inner tube as the 
Lennard- Jones potential ii(ry) = 4e[(cr / 'r^) 12 — (cr/r^) 6 ], where is the distance between 
atom i in the inner tube and atom j in the outer tube, e = 28 K, and a = 3.4AI19I. We 

n 

calculate the potential V(9) [Tj| when the outer tube rotates around the inner tube with 
angle 6 and plot it in Figdl We find that V{6) is periodic (with period 7r/2) and violates 



the spatial inversion symmetry. Thus the second condition to make a molecular motor is 
satisfied. Up to now we can explain why we select (8,4) tube and (14,8) tube. The following 
three criteria are considered: (1) layer distance is about 3.4A; (2) the shape of V{9) is not 
too complicated; (3) the difference of maximum and minimum of V(9) is remarkable. We 
find that only (8,4) and (14,8) tube satisfy these criteria for < ni,n 2 < 30 through our 
calculation. 

The easiest way to satisfy the first condition is to put our system into a thermal bath 
full of He gas whose temperature varies with time. We will show that the outer tube will 
exhibit a directional rotation. 

We can write the Langevin equation for the outer tube mp 2 9 = —V'(9) — r\9 + £(£), 
where m and p are respectively the mass and the radius of the outer tube, rj is the rotating 



viscous coefficient 2l(, and dot and prime indicate, respectively, differentiations with respect 
to tune ( and ang.e B. «t) is thenna. noise whieh satisfies ({(,)> = and the donation- 
dissipation relation (£(£)£(s)) = 2r]T(t)5(t — s) [22], where T(t) is temperature and the 
Boltzmann factor is set to 1. Let us consider the overdamped case that the inertial term 
mp 2 9 is much less than thermal fluctuations and can be neglected. Thus we arrive at 
T]9 = —V'(9) + and the corresponding Fokker-Planck equation [2I I23I] 



dP(9,t) _ d 
dt ~ 09 



V'(9)P(9,t) 



T(t) d 2 P(9,t) 
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where P(9, t) represents the probability of finding the outer tube in angle 9 at time t which 
satisfies P{9 + 7r/2, t) = P(9,t). If the period of temperature variation is T, we arrive at 
the average angular velocity in the long-time limit 

1 f t+T ,. r'\.\ V'(9)P(9,t) 



1 ft+-l r n/2 r 

(9) = lim - / dt d9\- 
t^oo T J t J [ 



(4) 

V J 

For example, set T(t) = T[l + A sm(27vt/T)} with T = 50 K and \A\ < 1. Let D = r}/f, 
t = Dt, U(9) = V(9)/T, T = DJ and P(9,t) = P(9,Dt), We arrive at the dimensionless 
equations of Eqs. (jHJ) and (jlj) 

BP 8 ~ 1-kt d 2 P 

^^ r(W + (1 + .4s m ^, (5, 

d9 1 r +J r /2 r ~i 

(^-) = lim - / dr d9 -U'{9)P . (6) 



x dr T-aoo J 

Let A = 0.01, we can numerically solve Eq.© and calculate Eq.(JH} with different J . The 
solid line in Fig0] shows the relation between the average dimensionless angular velocity (M-) 



and the dimensionless period J of temperature variation. We find that ~ for very 
small and large J y and ^ for middle J y which implies that the outer tube has a 
evident directional rotation in this period range. Thus we have constructed a temperature 
ratchet. 

For He gas at temperature T = 50 K, we can calculate rj = 861 Kns from its value at 
273 K 2lll3] and D = 17.2 ns. Based on these data, we obtain (jg) = -43 mad when 
J = 0.17, i.e. (9) = ^(jr) = —2.5 nrad/ns when T = 2.9 ns. Here the minus sign of average 
angular velocity means that the rotation of the outer tube around z-axis is left-handed. We 
notice that (8) = —2.5 nrad/ns is a remarkable value (-2.5 rad/s) that is easy to be observed 
from experiment. If consider the inertial effect of the outer tube, this value may be changed. 
But we believe that it is still an observable quantity. 

Furthermore, from Fig. Q]we observe that the sign of changes from " + " to " — " and 
back to " + " with J increasing, which suggests that the outer tube's rotation around the 
inner tube is the right-handed for small and large period J7, and the left-handed for middle 



period J . Through detailed analysis, we find that (^) is proportional to J' 3 for very small 
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period J and proportional to J for period J large enough, which agrees with the results 
of asymptotic analysis 

To summarize, we have theoretically constructed a temperature ratchet using a DWNT, 
but some practical difficulties are remained: How to synthesize the proper DWNT? How to 
forbid the axial translation of outer tube? If these obstacles are ruled out, molecular motor 
of DWNT must be reality. 

We are grateful to H. J. Zhou for reading the manuscript and giving much constructive 
advice. We thank to R. An for mathematical discussion. 
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FIG. 1: The unrolled honeycomb lattice of a SWNT. By rolling up the sheet such that the point Aq 
coincides with the point corresponding to chiral vector C^, a nanotube is formed. Vectors ai and 
a2 are real space unit vectors of the hexagonal lattice. The translational vector T is perpendicular 
to Ch and runs in the direction of tube axis. The vector R is the symmetry vector. Aq, Bq and 
A[, Bi{l = 1, 2, • • • , N) are used to denote the sites of carbon atoms. 




FIG. 2: A double-walled carbon nanotube with the inner tube's index being (8, 4) and the outer 
tube's index being (14, 8). z-axis is the tube axis, x-axis perpendicular to z passes through carbon 
atom Aq (see Fig^) of the inner tube and y-axis is perpendicular to the xz-plane. There is no 
obviously relative motion along radial direction between the inner and the outer tubes at low 
temperature. If we forbid the motion of the outer tube in the direction of z-axis, only the rotation 
of the outer tube around the inner tube is permitted. 
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FIG. 3: The potentials V(9) between the outer and the inner tubes with the outer tube rotating 
around the inner tube. 9 is the rotating angle. Here we have set V(0) = 0. The squares are the 
numerical results which can be well fitted by V{9) = 15.7 — 0.6 cos 49 H — 2. 2 sin 49 — 12.7cos 89 — 
6 sin 86>-1.7 cos 120 + 10.8 sin 126> (solid curve). 




FIG. 4: The average dimensionless angular velocity {d9/dr) of the outer tube rotating around the 
inner tube in thermal bath whose temperature changes with the dimensionless period J. The 
minus sign means the rotation around z-axis is the left-handed. 
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